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Abstract 

We consider the propagation and damping of isovector excitations 
in heated nuclear matter within the Landau Fermi-liquid theory. Re- 
sults obtained for nuclear matter are applied to calculate the Giant 
Dipole Resonance (GDR) at finite temperature in heavy spherical nu- 
clei within Steinwedel and Jensen model. 

The centroid energy of the GDR slightly decreases with increasing 
temperature and the width increases as T 2 for temperatures T < 5 
MeV in agreement with recent experimental data for GDR in 208 Pb 
and 120 Sn. 

The validity of the method for other Fermi fluids is finally sug- 
gested. 



PACS numbers: 24.30. Cz, 21.60.Ev, 21.65.+f 



In recent years the GDR built on highly excited states is in the center of 
many experimental and theoretical studies (c.f. |l| and references therein). 
In this context, one of the most important open problems is the behaviour 
of the GDR width in nonrotating nuclei as a function of temperature. There 
are two essentially different theoretical approaches to this problem. The first 
one || explains the temperature increasing of the width as an effect of the 
adiabatic coupling of the GDR to thermal shape deformations. In the second 
approach [Q, [|, |IJ the thermal contribution to the damping width arises from 
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an increasing nucleon-nucleon collision rate (2p2h excitations) plus a Landau 
spreading due to thermally allowed ph transitions || |7], || ||. 

In the present work, following the ideology of the second approach, we 
consider isovector volume vibrations in spin-isospin symmetrical nuclear mat- 
ter at finite temperature. A similar problem was considered in Refs. |7], § 
within the RPA method. However the Landau damping mechanism of the 
dissipation of a propagating mode due to thermal smearing of Fermi distribu- 
tion is too weak to be responsible for the fast increase of the observed GDR 
width with temperature [0, |[ |TTJ. This problem can be solved by taking 
into account the two-body dissipation through the collision integral of the 
Landau- Vlasov equation ||. The use of a quantum kinetic equation leads to 
memory effects in the collision term in order to include off energy-shell con- 
tributions ||12|| . Moreover it was shown in Refs. |13| , 14j] , that memory effects 



are essentially increasing the widths of multipole resonances at small tem- 
peratures. In this Rapid Communication we calculate the isovector strength 
function of nuclear matter taking into account both thermal Landau damping 
and two-body collisional dissipation, including the quantum memory contri- 
bution. 

The isovector response of uniform nuclear matter is described by the 
linearized Landau- Vlasov equation with a collision term treated in the relax- 



ation time approximation 0, 0, 14 



JU/ + v ■ V r 6f - V r (8U + 26V) ■ V p f cq = -~<5/|,>i , (1) 

where 5f = Sf n —6f p and 5U = 5U n —5U p are differences between neutron and 
proton distribution functions (d.f.) and mean fields respectively, 5V = 5V n — 
SV P is external field (^V^ = r q 5V, r n = +1, r p = —1) |7j, f eq (e p = p 2 /2m) is 
the equilibrium finite temperature Fermi distribution, and the notation / > 1 
means that the perturbation of the d.f. Sf\i>i in collision integral includes 
only Fermi surface distortions with multipolarity I > 1 in order to conserve 
the particle number in collision processes The inclusion of the I = 1 



harmonic in the collision integral of Eq. (|I|), at variance with the isoscalar 
case 0, is due to nonconservation of the isovector current, i.e. due to a 
collisional friction force between counterstreaming neutron and proton flows. 
The dynamical component of the isovector mean field 5U can be expressed 
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in terms of the isovector Landau parameter Fq : 

«r = CO 

where 

W'i«) = /^ i V(r,Pi«), (3) 
is the density perturbation, g = 2 is the spin degeneracy factor and 

^-/^(-^) w 

is the thermally averaged density of states, N(0) = gpFm/2ir 2 h 3 , where we 
put for simplicity m* = m = 938 MeV. 

For an external field 5V oc exp(ikr — iut), periodic in space and time, the 
isovector collective response function can be derived from Eq. dH): 



where Xt * s the intrinsic response function [|T|, [L6| . The explicit form of the 
function x^(u;,k) is (details of derivation in Ref. M): 

Xt{ ' ~ mp F N(T) J P s' + ts" X (ps/p) de p ' 1 ' 



where 



P = Pf[ — ), (7) 



5 f gdp 



3p eq 7 (27r^) 



e P /eq(e P ) , (8) 



are quasiparticle average momentum, average kinetic energy ( respectively 
normalized at T = on p F and e F ) and density, with the complex variable 



3 



s = s' + is", s" = — , s' = — , 10 
rpk pk 

x{z) is a Legendre function of the second kind 



-i 



Eq. (^D for the intrinsic response function of isovector vibrations has only a 
minor difference with isoscalar case. Namely, to recover the isoscalar response 
function given by Eq. (30) in Ref. 0, one should change is" — *> is"(l + 
3s'sp 2 /p 2 ) in the denominator of Eq. (^). This difference is just due to 
inclusion of the damping of the I = 1 harmonic in the isovector channel. 
We note, that the r.h.s. of Eq. (30) in Ref. || has an error: it should be 
multiplied by a minus sign. 

For a given momentum transfer k, the strength function per unit volume 

is: 

S t M--Im(x )~ {l + F , mx r T)? + {F!ilm{x r TW - (12) 

The strength function satisfies the following energy weighted sum rule 
(EWSR)il: 

7 k 2 

J duuS k {u) = — p , (13) 
o 

where po = 0.16 fm -3 is the nuclear saturation density. 

Collective modes are given by poles of the response function (|5|): 

1 + FUt(s)=0. (14) 
By solving Eq. (|H]) we obtain the complex frequency : 

uj — ujr + iui = k—(s — is") . (15) 
m 

The application of the formalism discussed above to finite nuclei is based 
on the Steinwedel- Jensen (SJ) picture JF], H [L7| which describes the GDR 



in heavy nuclei as a volume polarization mode conserving the total density 
P — Pn + p p - According to this model, we choose the wavenumber of the 
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normal mode as k = n/2R, where R is the radius of a nucleus. Inside the 
nucleus, the unperturbed distribution of nucleons is supposed to be uniform. 
The SJ model gives a good overall reproduction of the ground state GDR 
energies for heavy spherical nuclei 



We have to remark that the calculation of GDR widths in finite nuclei 
is a much more difficult problem. It can be only partially solved within 
our nuclear matter approach, since shell effects and the escape width due 
to particle emission are not taken into account in the present investigation. 
However in the temperature region T = 1 -v- 3 MeV the shell effects start 
to smear out and the escape width is still very small. Thus we expect our 
calculations to be quite reliable in this temperature region. 

Eq. ([!]) contains two free parameters: the isovector Landau parameter Fq 
and the two-body relaxation time r. 

The Landau parameter Fq at zero temperature can be expressed as a 
function of the symmetry energy coefficient (3 in the Weizsacker mass formula 
as follows jl9fl : 

F '(T = 0) = ^-1. (16) 

For the standard value of (3 = 28 MeV we have Fq(0) = 1.33. The value of 
the Fq decreases with temperature due to the decrease of the level density 
N(T): 

7T 2 /T\ 21 



F^T) ~ iftO) 



1 12 \e F ) 



(17) 



However the coupling constant Fq(T)/N(T) is independent on temperature. 

The two-body relaxation time r includes temperature and memory ef- 
fects : 

T = T 2 + ( W2tt) 2 ' ( } 

The dependence of the relaxation time on the frequency ujr arises from mem- 
ory effects and corresponds to the Landau prescription [[12]. The coefficient 



depends on nucleon-nucleon scattering cross sections. We have calcu- 
lated this coefficient using energy and angular dependent differential cross 
sections of pp and np scattering derived from Bonn A potential both with 



and without in-medium corrections |20], |21| (see Appendix). Results are : 
a* - - 1 = 2.3 (5.5) MeV in the case of vacuum (in-medium reduced) cross sec- 
tions. 
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Fig. 1 shows the photoabsorption cross section by a thermally excited 
nucleus 208 Pb, which can be expressed in terms of the strength function ([12]) 
as follows: 

_ . 47r 2 e 2 NZ + n . . 
a abs {huj) = ^-p- — — ftu;S fc (u;) . (19) 

Expression ( |i~9|) is obtained from comparison of the EWSR (|I3|) and the 
dipole sum rule [O 



/ 



o 



As the temperature is growing, the centroid energy of the GDR (i.e. the peak 
of the photoabsorption cross section) is slightly shifting to the left and the 
width is increasing. At variance with pure mean field predictions J7], ^ |], 
we see a shift to low frequencies due to the collision integral in Eq. (|]) (c.f. 



Refs. [£2], £3], Indeed at larger temperatures an increasing two-body 

dissipation should reduce the frequency of the collective motion, in close 
analogy with a classical oscillator with a friction force. 

In Fig. 2a (solid lines) we report the temperature dependence of the full 
width half maximum (FWHM) and the centroid energy Eqdr as functions 
of temperature, calculated from the photoabsorption cross section (|19"P for 
208 Pb. In parallel, in Fig. 3 (solid lines) the real and imaginary parts of the 
pole of the response function ([|) are shown. As far as the collective mode is 
underdamped, i.e. |Im(o;)|/Re(a;) -C 1, an approximate relation 

FWHM = 2|Im(w)| (21) 

has to be fulfilled jnj. We see from comparison of Figs. 2a and 3, that 
the condition fl2"Ip is really satisfied, corresponding to a Breit-Wigner-like 
shape of the photoabsorption strength (see Fig. 1). However, at higher 
temperatures T > 4 MeV the photoabsorption strength becomes much more 
closer to the Lorentzian shape p5| . 



At low temperatures, one can neglect the temperature spreading of the 
equilibrium Fermi distribution substituting 5(ep — e p ) instead of 
(—df eq (e p )/de p ) into Eq. (|). Thus we obtain the following approximate low 
temperature expression for the intrinsic response function : 

Xt(s) - T^#7T > (22) 
s' + is"x(s) 
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where variables s, s" and s' are defined by Eqs.flTOD with the change p — > 
Pf- We remark that Eq.(|^) when applied to the case of an electron gas 
(F = N(T) Ave 2 /k 2 ) gives just the longitudinal dielectric function 

e(k,w) = 1 + F x t t (23) 
obtained by Mermin |26|1 . This is to stress the more general framework of 



the results presented here [^7[ . 

In the rare collision regime (urt ^> 1), an approximate solution of the 
dispersion relation ( |i~4"|) with the intrinsic response function Eq.(p2|) can be 



found analytically (c.f. Refs. [12], £8] 



u R ~ v F ks {0) + O (T 4 ) , (24) 

^ (2F- + 1)((,W) 2 -1)-(F-) 2 

" 7 " rF&Fl - Wf + 1) ( )? ( } 

where s^ - 1 is the root of collisionless dispersion relation 

l+F> X (s^) = 0. (26) 

The Landau parameter Fq in Eqs.(25), fl26|) is taken at T = 0. The simple 



expression Eq.(23) for imaginary part of the frequency oo (long-dashed line 
on Fig. 3) reproduces the results of a numerical solution of the "exact" dis- 
persion relation (|TJ) (solid line on the same Fig. 3) with a good accuracy for 
temperatures T < 2 MeV. At larger temperatures, a slight increase of the 
damping due to temperature smearing of the Fermi distribution is obtained 
with the dispersion relation Eq. (|i~4"|) . The difference between these two so- 
lutions is of the order of the Landau damping rate in the pure mean field 
approach 0, [| . We see that thermal Landau damping is small for the case of 
GDR, at variance with the case of the breathing mode 0, since the isovec- 
tor Landau parameter is larger than the isoscalar one for nuclear effective 
interactions at normal density. As a consequence a weaker coupling between 
single particle and collective motion is expected for isovector vibrations |0J . 

The dominance of the collisional contribution to the total damping rate 
of the GDR is just expressed by an approximate relation (short-dashed line 
on Fig. 3) 

- Im(w) ~ - , (27) 
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which was used, for instance, in Ref. [13j to calculate widths of giant res- 
onances. The main deviation from formula ( p7|) in the dispersion relation 
Eq.(TH|) is caused by the exclusion of the / = harmonic from the collision 
integral in the r.h.s. of Eq. ([!]), i.e. due to taking into account the particle 
number conservation. This results in a smaller absolute value of the Im(u;) 
(see Fig. 3). 

As already discussed, a source of uncertainty in our calculations is given 
by the choice of the nucleon-nucleon cross sections. In Fig. 2 we present 
calculations with in-medium reduced cross sections (solid lines) and with free- 
space cross sections (dashed lines) in comparison with experimental widths 
for nuclei 208 Pb (a) and 120 Sn (b). There is a quite good reproduction of the 
experimental trend independently on the choice of cross sections. 

In conclusion, we have studied the isovector response of heated spin- 
isospin symmetric nuclear matter on the basis of the linearized Landau- 
Vlasov equation with a collision integral including memory effects in the re- 
laxation time approximation. The contribution of the thermal Landau damp- 
ing into the total damping rate of isovector vibrations is found to be very 
small. Thus the relaxation of the volume isovector mode is caused mainly by 
nucleon-nucleon collisions. Increasing temperature shifts the centroid energy 
of the isovector strength function to smaller values. The calculated width is 
proportional to 1/r in the temperature region T < 5 MeV studied in this 
letter, that corresponds to the collisional damping of the isovector zero sound 
mode. This leads to a T 2 behaviour of the GDR — FWHM in the region 
where the collective mode can propagate. 

We have shown that some general GDR properties at high excitation en- 
ergy can be obtained directly from Fermi liquid theory. However the aim of 
this work is not to get a perfect agreement with data for finite nuclei particu- 
larly at low temperatures. Indeed we are well aware that other contributions 
to the damping are missing in the present approach: (i) The fragmentation 
width, observed in RPA calculations, which can be interpreted as a Landau 
damping mechanism in finite Fermi systems present also at zero tempera- 
ture (ii) Thermal shape fluctuations |2|, (iii) Fluctuations due to 
nucleon-nucleon correlations |33] . 



A recent statistical model analysis of 7-spectra produced by inelastic 
a-scattering on 120 Sn p5f resulted in conclusion that neither the thermal 



fluctuation model of Ref. |39[ nor the collisional damping model of this work 
could reproduce data in details. We believe that the combination of the two 
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models would give a much better agreement. 

Extension to isospin asymmetric nuclear systems, more suitable for the 
Pb case, can be performed following the approach of ref.s |B3, B3, We 



do not expect to have substantial variations in the temperature behaviour of 
the isovector and isoscalar modes unless very high charge asymmetries are 
reached |33, 35 ] . However in charge asymmetric nuclei, new soft mode, differ- 
ent from isovector and isoscalar ones, appears [EH] due to collisional coupling 



of proton and neutron vibrations, that requires further investigations. 

Finally a comment on the validity of the formalism applied here to the 
study of hot GDR in nuclei to a broader context of physical problems. In- 
deed it can be easily generalized to the description of relative vibrations of 
any two-component Fermi liquid with a mutual attraction, for instance of a 
Coulomb plasma consisting of opposite charged fermions. Another applica- 
tion could be to the oscillations of the electronic cloud in metallic clusters, 
where momentum nonconserving I = 1 term in collision integral appears due 
to scattering of electrons on impurities [p6fl. 
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Appendix 

Here we will calculate the coefficient in Eq. (pl|) for the collisional 
relaxation time r. For simplicity, we will use in derivation the Boltzmann- 
Uehling-Uehlenbeck (BUU) collision integrals without memory effects. How- 
ever, as it is shown in Ref. (|^6|), the collision integrals with memory 
effects give the result which can be obtained using Landau prescription 
t — t buu/(1 + (^ a; i?/27rT) 2 ), where tbuu — ha^/T 2 is the relaxation time 
given by BUU collision integrals. 
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Time evolution of the space-uniform isovector d.f. / = f n — f p satisfies 
the equation (gj| |14|) : 

dfiP't) T J ,T _J _J /OON 

1 1 nn i 1 np 1 pp 1 pn > \ °) 

where I qiq2 stands for the collision integral of particles of the sort q\ = n,p 
with particles of the sort q 2 = n, p. Explicitly : 

4 f 

J <?i? 2 (Pi>*) = ( 2 7T^) 6 J rf P 2rf P 3rf P 4U, <?i92(Pi>P2;P3,P4)5(Ae)5(Ap) 

Q(U(Px), /»(Pa); /«(Pa), / ?2 (P4)) , (29) 

where 

Wi, / 2 ; /a, / 4 ) = (i - /i)(i - h)hh - A/ a (i - / 8 )(i - / 4 ) ; 

u 'gig2(Pi) P2; P3? P 4 ) is the spin- averaged probability of two-body collisions 
with initial momenta (pi,P2) and final momenta (p3,P4); Ae = t\ + e 2 — 
£3 - e 4 , e i = Pi /2m, i = 1, 2, 3, 4, Ap = p x + p 2 - p 3 - p 4 . Neglecting the 
dependence of w qiq2 on the d.f., we can write down perturbations of collision 
integrals (|29|) keeping terms of the first order in 8f q : 

4 /■ 

^glga(Pl^) = ( 2?r ^)6 7 rf P2^P3rfp 4 W gi(?2 (pi,P2;P3,P 4 )5(Ap) 

{a (1) Vffi(Pi) + « (2) ^ 2 ( P2 ) W 3) W(Ps) +« (4) ^ 2 (P 4 )} , (30) 

where 

MP) = ^ ? (P) (^) , (ff = n,p), (31) 
(0 = 5Q(feg(ei), feg(e 2 ); feq(e 3 ), fegfa)) 

&feq(ei) 



dei 



-S(Ae) , 2 = 1,2,3,4. (32) 



In Eq. (^) the isospin-symmetric nuclear matter is considered that results 
in the same equilibrium d.f. for neutrons and protons. 
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For the perturbation of the collision integral I of the r.h.s. of Eq. (|28|). 
assuming isotopic invariance {w pp = w nn , w pn = w np ) we can write after 
simple algebra : 

81 5I nn -\- 8I np 8I pp 8Ipn 
4 r 

I dp 2 dp 3 dp 4 5(Ap){ (w pp + w np ) V(Pi) + " (3) V'(Pa)) + 
W „ p )(« (2) ^(p 2 ) + « (4 V(P4)) } , (33) 



where ijj(p) = ip n (p) ~ ^p(p)- 

The triple integral over momenta in (RBf) can be taken using Abrikosov- 



Khalatnikov transformation ||T2j| , which is valid in the limit T -C ep ■ 

r (tyi ) r sin fc) r f f 
J dp 2 dp 3 dp 4 S(Ap) = — - — J dQ / d*Jd<h de 2 de 3 de 4 , 

^ ' ' 

_ (34) 
where 9 = (pi,p 2 ) is the angle between momenta of colliding particles (p; = 
Pi/|Pi|) i = 1)2,3,4); is the angle between planes given by momenta of 
incoming and outcoming particles : 

[pi x p 2 ] • [P3 x pi] . . 

cos0 = — — — — ; 35 

Ibi XP2]||[p 3 Xp 4 ]| 

02 is azimutal angle of p 2 in the system with z-axis along pi. 

We decompose the perturbation of the d.f. into spherical harmonics : 

ip(pi,t) = ®im(Pi,t)Yi m (pi) , i = 1,2,3,4, (36) 

Ijn 

where coefficients ai m can be taken on the Fermi surface, since T <C ep. 

According to Refs. fi3| , |37 |, the partial relaxation time t\ is defined 
as follows : 

oo 

1 o 



Bl/t/ 



r, 



/de/dty^Cp) Sf(p,t) 



JdQ p Y l * m (p)5I(p,t) ^ 
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where 



SI(p) = / deSI(p) , e = p 2 /2m 



(38) 



Using ([34"1), after somewhat lengthy but standard calculations, we come 
to the expression : 



1 



' ' {< W pp ^ +) > +2 < Wnp&i ] >} = — 
*• - 1 Ki 



(39) 



where angular brackets denote the averaging over angles 9 and [T||, [28 

/ ,IH — 

T J cos „ 

2 



<F(9,<j>) >=^- fdB^ ld<pF{6 

lit J COS 77 •/ 



(40) 



(+) 
i 

(-) 



1 + i^(p2 ■ Pi) - ^l(P3 ■ Pi) - PliPi ■ Pi) , 
1 - Pl(j>2 ■ Pi) - fl(Pa • Pi) + Pl{p4, ■ Pi) • 



(41) 

(42) 



A factor 2 at the second term in curly brackets of (|3TS|) is due to half mo- 



mentum space integration over <ip 3 in the l.h.s. of Eq. (|34|) |12| , p8| . The 
arguments of the Legendre polynomials in functions (|4"T|) , (f4"2"|) are : 



P2 - Pi 
P3 "Pi 

P4 ■ Pi 

For I = 1 , 2 and oo we have 



cos 8 , 

2 . o 9 
cos — h sin - cos < 
2 2 

cos sin - cos ( 

2 2 



= 



(-) 



(+) 



3 sin - sin 
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2 " = 3 Slit* 

$(+) = $(-) = ! _ 



W • 2 17 • 2 V 

4 sin -sin — , 

2 2 ' 

$jf } = 3 sin 2 9 sin 2 ^ , 



(43) 

(44) 
(45) 
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Collision probabilities can be expressed in terms of cross sections as fol- 
lows : 



Wpp = (2^^L ; (46) 
Wpn = (2tt^) 3 da pn ^ (47) 

where da pp / dfl c , m , and da pn / dfl c , m , are differential cross sections of pp and 
np scattering; dVL cm ^ = sin9 c . m .d6 crn ,d(f) c . m /, Q c . m . and c . m . are polar and 
azimutal scattering angles in the center of mass system of colliding particles; 
/i = m/2 is the reduced mass. Differential cross sections da pp / dfl c , m , and 
d(J pn / dfl c , m , depend on the relative momentum pf = \p 1 — P2I/2 of scattered 
particles and on the polar angle 

a Api -P2) • (Pa -P4)\ 
v r m — arccos — - — 



V iPl - P2||P3 - P4| / 

For particles scattered on the Fermi surface, we have : 

Q 

pf = p F sin - , 9 c . m . = . 

In a particular case of isotropic energy independent cross sections the result 
of Ref. is recovered : 

1 32ma„ T2) (4g) 



Tl BUU " 9 h 3 



32m °°T\ (49) 



t? uu 15 h 

where a v = a np /2, a s = (a nn + cr pp + 2a np )/4, a np = (Aii)da np / dVt c , m , ~ 50 
mb, a nn ~ a pp = {2-n)da pp / dVt c , m , ~ 25 mb. 

We derived relaxation times t buv , r^ c/c/ and t^ uu using pp and nn en- 
ergy and angular dependent differential cross sections calculated with Bonn 



A potential |20], [21|| . Results of these calculations both with vacuum and 
in-medium cross sections at normal nuclear density are given in the Table. 
It is seen from the Table that always Ti UU ~ T2 UU — t^ uu . This gives an 
idea to put the same value for all relaxation times t® vij ' , i.e. to apply usual 
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relaxation time approximation. Thus, we obtain the collision integral of Eq. 
([!]) with relaxation time r given by Eq. (18), where 



a 



(-) 



h 



2.3 (5.4) MeV 



for vacuum (in-medium) cross sections of Ref. pT |. 
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Table Parameters kj, 1 — 1,2 and oo (MeV 2 fm/c) defined in Eq. ( p9[ ) at 
various choises of nucleon-nucleon scattering cross sections. 



NN cross sections 




k 2 




Vacuum of Kef. 12 1| 


503 


491 


401 


In-medium at 

P = Po 
ofRef. £1 


1123 


1068 


1003 


Isotropic energy- 
independent 
ofRef. || 


920 


1022 


818 
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Figure captions 



Fig. 1 Photoabsorption cross section by an excited nucleus 208 Pb as a func- 
tion of photon energy for temperatures T = MeV (solid line), T = 2 
MeV (short-dashed line) and T = 4 MeV (long-dashed line) calculated 
with in-medium cross sections of Ref . . 

Fig. 2 Centroid energy E GDR (MeV) and FWHM (MeV) for the GDR 
mode in nuclei 208 Pb (a) and 120 Sn (b) as functions of temperature 
calculated with in-medium cross sections (solid lines) and with free 
space cross sections (dashed lines). Points with errorbars show the 
experimental widths from Ref. P3 - 208 Pb and from Ref. |30] - 120 Sn. 



Fig. 3 Temperature dependence of real and imaginary parts of the pole of 
the response function (||D - upper and lower solid lines respectively, and 
collisional width T co u = 2h/r - short-dashed line. The long-dashed line 
shows the imaginary part of pole as given by the approximate Eq . (|25|) . 
All values are in MeV. In-medium cross sections are used. 
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